The electron transport through a quantum dot in the Coulomb blockade regime: 
Non-equilibrium Green's functions based model 
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We explore electron transport through a quantum dot coupled to the source and drain charge 
reservoirs. We trace the transition from the Coulomb blockade regime to Kondo regime in the 
electron transport through the dot occuring when we gradually strengthen the coupling of the 
dot to the charge reservoirs. The current-voltage (/ — V) characteristics are calculated using 
the equations of motion approach within the nonequilibrium Green's functions formalism (NEGF) 
beyond the Hartree-Fock approximation. We show that within the Coulomb blockade regime the 
I — V characteristics for the quantum dot containing a single spin-degenerated level with the energy 
Eo include two steps whose locations are determined by the values of Eo and the energy of Coulomb 
interaction of electrons in the dot U. The heights of the steps are related as 2 : 1 which is consistent 
with the results obtained by means of the transition rate equations. 

PACS numbers: 73.23.Hk, 72.10.Di, 73.63.-b 



I. INTRODUCTION: 



As known, quantum dots host such interesting quan- 
tum transport phenomena as Coulomb blockade [l| and 
Kondo effect 0, Fast development of experimen- 
tal technique in recent years provides broad opportu- 
nities for fabrication of devices including novel active 
elements such as single molecules, nanowires and car- 
bon nanotubcs, placed in between metallic leads. Such 
junctions are commonly treated as quantum dots cou- 
pled to the macroscopic charge reservoirs which play 
the parts of source and drain in the electron transport. 
Current achievements in nanoelectronics make it possi- 
ble to observe electron transport through realistic junc- 
tions varying external parameters (source-drain and gate 
voltage and temperature) and internal characteristics of 
the junction, e.g. the coupling strength of the quan- 
tum dot to the source and drain reservoirs. As a re- 
sult, both supression of the electron transport through 
the quantum dot at small values of the source-drain volt- 
age (V) arising due to the high charging energy in- 
side the dot (Coulomb blockade), and the increase in 
the electrical conduction of the junction at V w due 
to the Kondo effect were repeatedly observed in exper- 
iments on molecular and carbon nanotubes junctions 

The Coulomb blockade occurs when the quantum dot 
representing a molecule or a carbon nanotube is weakly 
coupled to the source and drain which means that the 
charging energies U are much greater than the contact 
and temperature induced broadening of the electron lev- 
els in the dot. Theoretical studies of the electron trans- 
port through a quantum dot within the Coulomb block- 
ade regime mostly employ "master" equations for the 
transition rates between the states of the dot differing by 



a single electron 15, 3, 17, HI]. It was shown that the 
current-voltage curves for a dot including a single spin- 
degenerated level exhibite two steps. The steps corre- 
spond to the gradual adding/removing two electrons with 
the different spin orientations to/from the dot. Assuming 
the symmetric contacts of the dot with both charge reser- 
voirs, the first step in the current-voltage curve must be 
two times higher than the second one (seeRefs. [II [H]). 
This result has a clear physical sense. Obviously, one may 
put the first electron of an arbitrary spin orientation to 
the empty level in the dot. However, the spin of the 
second electron added to the dot is already determined 
by the spin of the first one. So, there are two ways to 
add/remove the first electron to/from the empty/filled 
level in the dot but only one way to add/remove the sec- 
ond one. 

While the approach based on the transition rates "mas- 
ter" equations brings sound results within the Coloumb 
blockade regime, its generalization to the case of stronger 
coupled junctions is not justified. More sophisticated 
methods employing various modifications of the nonequi- 
librium Green's functions formalism (NEGF) are devel- 
oped to analyze electron transport through quantum 
dots coupled to the leads [13] ■ These theoretical tech- 
niques were repeatedly employed to calculate the re- 
tarded/advanced and lesser Green's functions for the dot. 
The Green's functions are needed to compute the elec- 
tron density of states (DOS) on the dot where the Kondo 
peak could be manifested, as well as the electron current 
through the junction. A commonly used approach within 
the NEGF is the equations of motion (EOM) method 
which was first proposed to describe quantum dot sys- 
tems by Meir, Wingreen and Lee [2l|. This technique 
enables to express the relevant Green's functions in terms 
of higher-order Green's functions. Writing out EOM for 
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the latter, one arrives at the infinite sequence of the equa- 
tions successively involving Green's functions of higher 
orders. 

To get expressions for the necessary Green's functions 
this system of EOM must be truncated. Also, higher- 
order Green's functions still included in the remaining 
EOM must be approximated to express them by means 
of the lower-order Green's functions. In outline, this pro- 
cedure is well-known and commonly used. However, in 
numerous existing papers important details of the above 
procedure vary which brings different approximations for 
the relevant Green's functions. Accordingly, the results 
of theoretical studies of the quantum dot response vary. 

The NEGF formalism was successfully applied to de- 
scribe the quantum dot response within the Coulomb 
blockade regime (weak coupling of the dot to the electron 
reservoirs), using the expressions for both retarded and 
lesser Green's function within the Hartree-Fock approx- 
imation [z| . The results for the electron transport 
through the quantum dot obtained in these works give 
the correct ratio of the subsequent current steps, namely: 
2:1, assuming the dot to be symmetrically coupled to 
the leads. It is a common knowledge that within the 
Hartree-Fock approximation one cannot catch the Kondo 
peak. To quantitatively describe this effect one may need 
to include hybridization up to very high orders in calcula- 
tion of the lead-dot coupling parameters, and the studies 
aiming at better theoretical analysis of the Kondo effect 
started by Meir et al 24[ are still in progress (see. e.g. 
Refs. [USUI). 

Obviously, the Green's functions suitable to reveal the 
Kondo peak in the electron DOS on the dot must pro- 
vide the proper description of the transport through the 
dot within the Coulomb blockade regime when leads- 
dot couplings are weak. However, usually NEGF based 
results beyond the Hartree-Fock approximation success- 
fully employed to describe the Kondo effect and related 
phenomena, were not applied to the limiting case of 
the Coulomb blockade. In those works were such ap- 
plication was carried out, the correct ratio of heights 
of the subsequent current steps, namely 2 : 1, was 
not obtained. For instance, Muralidharan, Ghosh and 
Datta [HI reported the results on NEGF calculations 
which give equal heights of the two steps on the current- 
voltage curve within the Coulomb blockade regime, and 
Galperin, Nitzan and Ratner studies [2(| result in the 
heights ratio about 1.6 : 1. So, up to present there ex- 
ists a discrepancy between NEGF and master equations 
based results. This discrepancy may not be easily disre- 
garded for its existence makes questionable general and 
useful results based on the advanced techniques within 
the NEGF formalism. Recently, the above disagreement 
was discussed in the Ref. [jjj . 

The purpose of the present work is to clarify the is- 
sue, and to erase the above disconsistency. We apply 
equation of motion based approach to NEGF method to 



get a description of both Coulomb blockade and Kondo 
regimes in the electron transport through a quantum dot 
within the same computational procedure for the relevant 
Green's functions. We show that the assumed approxi- 
mations give correct ratio of heights of the steps in the 
current- voltage curves within the Coulomb blockade limit 
and allow to get the Kondo peak, as well. For simplicity, 
we omit from consideration dissipative effects, and we 
concentrate on a non-magnetic quantum dot coupled to 
non-ferromagnetic leads. Obtained results are valid for 
an arbitrary value of the Coulomb interaction on the dot, 
assuming that the latter is stronger than dot-leads cou- 
pling. Also, the analysis may be generalized to magnetic 
systems and to dots containing many levels provided that 
the charging energy is small compared to the adjacent 
levels separations. We remark that the present analysis 
is confined with using simplest possible approximations 
for the necessary Green's functions. Therefore we do not 
take into account corrections proposed in Ref. [26| to 
better approximate the Kondo peak. 



II. MAIN EQUATIONS 

We write the Hamiltonian of the junction including the 
dot and source and drain reservoirs (leads) as: 



H = Ht 



Hr 



Hr 



Hn 



(1) 



Here, the terms Hp (/? = L, R) corrrespond to the left 
and right leads and describe them within the noninter- 
acting particles approximation 



(2) 



where e r p a are the single-electron energies in the elec- 
trode (3 for the electron states r, a (the index a la- 
beling spin up and spin down electrons) , and ci^ , c r p a 
denote the creation and annihilation operators for the 
electrodes. We assume that energy levels e r p a in the 
electron reservoirs are uniformly spaced over the ranges 
corresponding to the electron conduction bands for both 
spin orientations. The spacing between two adjacent lev- 
els Sep a is supposed to be much smaller than the band- 
width Ap a (Sep*/ bp* ~ lCr 4 ). 

The term Hp describes the dot itself. Taking into 
account the electron-electron interaction, this term has 
the form: 



H D = E a c+c a + Uc^c^c^ci 



(3) 



where , c a are the creation and annihilation operators 
for the electrons in the dot; E a = Eq is the energy 
of the single dot level, the energy U characterizes the 
Coulomb interaction of electrons in the dot. We carry 
out calculations within the low temperature regime, so 
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U is supposed to significantly exceed the thermal energy 
kT. The last term in the Hamiltonian Eq. Q] describes 
the coupling of the reservoirs to the dot: 



r/3cr 



(4) 



where T r $ a are the coupling parameters describing the 
coupling of the r, a electron states belonging to the elec- 
trode (3 to the quantum dot. 

Now, we employ the equations of motion method 
to compute the retarded Green's function for the dot. 
Since no spin-flip processes on the dot are taken into 
account, and we assume the electron transport to be 
spin-conserving, we may separately introduce the Green's 
functions for each spin channel. The computational pro- 
cedure is described in the Appendix. We obtain: 



G?(E) = 



E-E -Z% 2 ~U(l-(n_ a )) 
(E-E Q - £ a)(£ -Eq-U- ££ 2 ) + UH l0 



(5) 

Here, (n-a) = (c^i^C—J), and self energy corrections 
are given by: 



T r f)o\ 



y 

^ E- e rl3a + %r\ ■ 



(6) 



£ lCT =5> r/3 ,_ CT | 2 /t 



r/3 



E-2E -U + e rP .- a + it] 



E - e r p.-a + if] . 

S 2cr = Vrfi-af I ^ ^ i""^ 

^ [E- 2E - U + e r p- a + it] 

1 



E - e r /3,- a + if] 

^02 = ^Oo- + ^2(7. 



(7) 



(8) 
(9) 



Here, ff a is the Fermi distribution function for the en- 
ergy e r p a and the chemical potential ]iP and 7] is a 
positive infinitesimal parameter. The expression[5]for the 
Green's function G r J(E) was first derived by Meir and 
Wingreen 21| . Later, the same expression was obtained 
in Ref. [25 1 and (assuming that the electron energy lev- 
els on the leads are spin degenerate) in the book [l9( . It 
was repeatedly used to qualitatively describe the Kondo 
peak in the electron DOS (see e.g. Refs. Esj). As 
shown in these works, the Kondo peak originates from 
the last term in the denominator of the Eq. [5[ At low 
temperatures the real part of the self-energy term 
diverges giving rise to the peak at /i L = fi R = E. Better 
approximations for the Green's function G r J(E) such as 
that derived in the Ref. [2(| may give better results for 
the shape and height of the Kondo peak. However, the 
very emergence of the latter is accounted for in the Eq. 

13 



The occupation numbers \ n aj = ( c i c <r) an d 
(n-a) = (dt^C—a) may be found solving the integral 
equations of the form: 



(n a ) 



1 

2^ 



dEImiG* 



(10) 



where the lesser Green's fucntion G< is introduced. The 
function G< obeys the Keldysh equation: 



G<(E)=G r J(E)X<G?(E) 



(11) 



where G® a (E) is the Fourier component of the advanced 
Green's function. The Fourier transform of the lesser 
Green's function G< is related to the Fourier transforms 
of the retarded and advanced Green's functions of the 
form (5[ so the quantities (?v) and (ri— CT ) appear in the 
integrands in the right hand sides of the equations [TOl for 
both spin directions. In further calculations we use the 
following expression for : 



(12) 



Consequently, we have 

G< = iG r J{E)GT\E) [T L a ti(E) + T*f?(E)] (13) 

where f^ R (E) are the Fermi distribution functions for 
the left /right reservoirs, and 



where 



^Oct — 



-21m S 



E 



i/3 

J 0fT 



E 



£r/3<r + if] 



(14) 



(15) 



Employing the expression [12] for E< we assume that 
the Coulomb interaction of the electrons on the dot does 
not affect the coupling of the latter to the leads, as well 
as it occurs within the Hartree-Fock approximation [22l |. 
Such assumption agrees with the way of estimating av- 
erages in the process of decoupling of high-order Green's 
functions which was used to derive Eq. [SJ It seems rea- 
sonable while the dot coupling to the leads is weaker 
than the characteristic energy of the Coulomb interac- 
tions on the dot. The proper choice of approximation 
for the lesser Green's function is very important for it 
directly affects the values of occupation numbers 
and, consequently, the heights of the peaks in the electron 
DOS in the Coulomb blockade regime. Also, the relative 
heights of the steps in the current-voltage characteris- 
tics depend on the occupation numbers. The employed 
form for G< (E) differs form those used in some previous 
works (see e.g. Refs. [20, UM) where extra terms arising 
due to the Coulomb interaction on the dot are inserted in 
the expression for £<. Certainly, one may expect such 
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terms to appear as corrections to the main approximation 
1121 However, these terms must become insignificant and 
negligible in the Coulomb blockade regime when the dot 
is weakly coupled to the leads. The expression for E< 
reported in Ref. [2o| does not satisfy this requirement. 

Substituting Eq. [5] into Eq. [TBI and inserting the result 
into Eq. [lO]we obtain the system of equations to find the 
mean occupation numbers of electrons on the dot: 

(n ±a ) = P± a + U(n Ta )Q± a + U 2 (n Ta ) 2 R ±a . (16) 
Here, 

P ±° = Y E / ?lj P AE)A ±a A+ a dE, (17) 

Q±a = y £ / TljP(E) (A ±a B+ a + A+ a B ±a ) dE, 

(18) 



^ = ^E/ rLtf(E)B±*B± a dE, 



(19) 



and the expressions for A a and B a are obtained using 
the Green's function El namely: 



A n 



E-E Q -U-H° G 



(E-E - Z Qa )(E -Eq-U- Eg- 2 ) + C/E lCT ' 



B„ = 



1 



(E-E - E 0CT )(S -Eq-U- Eg 2 ) + f/E 1(T 



(20) 



When considering a nonmagnetic system, the energy 
levels in the source and drain (and in the quantum dot, 
as well) are spin degenerated, so A a = A- a , B„ — 
B- a . Correspondingly the occupation numbers for both 
spin orientations coincide with each other. In this case 
( n cr) = ( n -a) are described with the equation 



(n CT ) 



1-UQ- - UQ) 2 - 4U 2 RP 



2U 2 R 



(21) 



When the electron correlations on the quantum dot are 
weak (E7 — ► 0), the result [21] is reduced to (n CT ) = P a - 

The Eq. [2T] is an important result for it gives an ex- 
plicit analytical expression for the level populations, so 
we do not need a time-consuming self-consistent itera- 
tive procedure to find the occupation numbers 27J. Now, 
we write the expression for the electric current I flow- 
ing through the junction when we apply the voltage V 
across the latter. This expression was derived by Jauho, 
Wingreen and Meir (see refs. [28|, [29j), and it has the 
form: 



iey [dE . L z 



If/*) (G r J - GT) 



(If - If) G<] 



Also, we could employ the equivalent expression: 

a 

-i {Y L J L a + If /f ) G r JG a a a (T L a - r*)] . (23) 

Further we assume the symmetrical voltage division: 
f^L,R = Ef ± \ e Vi an( i we P u t -Ef = 0. Also, in further 
calculations we introduce a gate potential V g , so that 



E o (V g ) = E (V g =0)+eV g 



(24) 



(22) 



The application of the gate voltage shifts the energy level 
in the dot, so the energies E and E + U become asym- 
metrically arranged with respect to E = Ep . This asym- 
metry is necessary to obtain two steps in the I —V curve 
assuming the symmetric voltage division. Otherwise, the 
contributions to the current from the vicinities of E = Eo 
and E = Eq + U become superposed, and only one step 
emerges. 



III. DISCUSSION 

To simplify further calculations we assume all coupling 
strengths to take on the same value, namely: r r p a = r. 
Now, we calculate the current through the junction in 
the limit of the weak coupling of the dot to the reser- 
voirs (r <C U). To carry out the calculations we need 
to know the nonequilibrium occupation numbers and we 
compute them using Eq. 1211 The occupation numbers 
are sensitive to the value of the applied voltage at small 
voltages as shown in the Fig. 1. At low values of the 
bias voltage (V < OAV) the assumed dot energy level 
(Eo — — 0.2eV) is situated below the chemical poten- 
tials for both leads. Therefore the dot is able to receive 
an electron but unable to transfer it to another lead. 
Accordingly, the average occupation number is close to 
unity and the electric current through the junction takes 
on values close to zero (see Fig. 2). At V — OAV the dot 
energy Eq crosses the chemical potential fi R , and the 
dot becomes active in electronic transport. Now, the elec- 
tron which arrive at the dot from one reservoir may leave 
it for another reservoir. This results in a pronounced de- 
crease in the average occupation on the dot accompanied 
by an increase in the current. One more change in both 
average electron occupation on the dot and the current 
through the junction occurs at V = 0.6V when the en- 
ergy Eq + U (£7 = 0.56^) crosses fj, R . At higher voltage 
all curves presented in the figures 1,2 level off. The cur- 
rent reachs its maximum value and the average number of 
electrons in the dot reachs its minimum. The minimum 
occupation number is noticeably less than one but its 
value is nonzero for electrons unceasingly travel through 
the junction. 

The current-voltage curves in the Fig. 2 show typi- 
cal Coulomb blockade features, namely, two steps whose 
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FIG. 1: Average occupation numbers in the quantum dot as 
function of the voltage applied across the junction. The curves 
are plotted assuming U = 0.5eV, kT = 0.00026eV, E (V g ) = 
-Q.2eV, t = OMeV (solid line), 0.02eV (dash-dotted line) 
and O.OieV (dashed line). 



heights are related as 2:1. So, we see that the results 
concerning the electron transport through a quantum dot 
obtained employing the transition rates equations may be 
quantitatively reproduced within the NEGF formalism 
when the retarded and lesser Green's functions are ap- 
proximated by Eqs. [SJ H31 and the explicit expression I2T1 
for the occupation numbers is used. Therefore, the dis- 
agreement discussed in the begining of the present work 
may be successfully erased, and the consistency between 
the NEGF formalism and the transitions rates equations 
in the description of the Coulomb blockade regime could 
be restored beyond the Hartree-Fock approximation. 

It is worthwhile to remark that particular values of the 
average occupation numbers on the dot are very respon- 
sive to the gate voltage value V g (the latter determines 
how the dot energy level is situated with respect to the 
Fermi levels of the leads in the absence of the voltage 
applied across them) and to the Coulomb interaction en- 
ergy U. Therefore, different values chosen for V g and U 
lead to different average occupation numbers. However, 
under various assumptions for the V g and U values, one 
quantity does not vary provided the symmetric coupling 
of the dot to the leads. This quantity is the relative 
height of the subsequent steps in the average occupation 
numbers of the electrons in the dot which are revealed as 
the voltage across the leads increases. As shown in the 
Fig. 1 this ratio is 2:1, exactly the same as for the 
subsequent current steps in the I — V curves in the Fig. 
2. And it is this ratio which ensures the correct shape of 
I — V curves. For instance, comparing figure 1 with the 



FIG. 2: Current through a junction within the Coulomb 
blockade regime. The curves are plotted assuming U = 
0.5eV, kT = 0.00026eV, E (V g ) = -0.2eV, r = 0.04eC, r = 
O.OleV (solid line), 0.02eC (dash-dotted line) and O.OAeV 
(dashed line). The factor Tq/t 2 is introduced to bring all 
I — V curves to the same scale. The ratio of heights of the 
two steps revealed in the curves equals 2:1. 

corresponding result reported on the work [22j |, one may 
see that the values of the occupation numbers given in 
Ref. [13] considerably differ from those obtained in the 
present work. Nevertheless, the ratio of the subsequent 
steps heights in the occupation numbers versus voltage 
curves is 2:1, and this provides for the same ratio of 
heights of the subsequent current steps. 

Also, the Eqs. El Q2J [19] could be applied to analyze 
the Kondo effect which is manifested at stronger coupling 
strengths of the dot to the charge reservoirs. The electron 
density of states (DOS) in the dot is given as 

D(E)=-2j2^[G r J(E)]. (25) 

a 

Using G r J(E) in the form [5] we may compute the elec- 
tron DOS. The results are presented in the Fig. 3 where 
the equilibrium DOS is shown for three values of the cou- 
pling strength r. For a sufficiently strong coupling of the 
dot to the source and drain reservoirs (r = 0.2eV) the 
sharp Kondo peak appears at E = 0, and the peaks at 
E = E and E = E + U are damped. At weaker cou- 
pling strength (r = O.leV, 0.07eV) the Kondo peak is 
reduced to a tiny feature but the maxima at E = Eq and 
E — Eq + U which determine the conductance within 
the Coulomb blockade regime emerge. The weaker is the 
coupling the higher become these peaks. 

The heights of the peaks differ. Technically, this dis- 
tinction originates from the fact that the value of m_ a j 



6 




FIG. 3: The equilibrium (V = 0) electron density of states 
in the quantum dot. The curves are plotted for r = 0.2eV 
(solid line) O.leV (dashed line) and O.OGeV (dash-dotted 
line) . The remaining parameters are the same as those in the 
figure 1. 



FIG. 4: Splitting of the Kondo peak in the electron density of 
states with increasing bias. The curves are plotted for eV = 
(solid line); 0.02eV (dashed line) , Q.QAeV (dash-dotted line) 
and O.OGeV (dotted line). The coupling strength r = 0.2eV. 
The remaining parameters are taken as in the figures 1,2. 



in the expression for the retarded Green's function (5) 
differs from 0.5. As shown in the book [l9[, assuming 
(ri-cr) = ( r v) = 0.5 one would get the Coulomb block- 
ade peaks of equal height but such assumption leads to 
the wrong result for the current. Namely, one would 
obtain two steps of equal height on the current voltage 
curves. The physical reason of this difference in the peak 
heights is the same as for the difference in heights of the 
steps on the current-voltage characteristics within the 
Coulomb blockade regime. The latter was discussed in 
the Introduction. When the voltage is applied across the 
junction the Kondo peak splits in two maxima whose 
heights are significantly smaller than the height of the 
original equilibrium Kondo peak and the greater is the 
voltage the lower these maxima become. This is shown in 
the Fig. 4. So, the present formalism sufficiently repro- 
duces the results of earlier works concerning the Kondo 
effect (see e.g. Refs. [HHH). 

Finally, in the present work we theoretically analyzed 
the electron transport through a single quantum dot cou- 
pled to the source and drain charge reservoirs. The anal- 
ysis was based on the NEGF formalism. The expression 
(5) for the retarded Green's function was obtained us- 
ing the equations of motion method, an d ag rees with the 
results of the previous works 2lj, [24|, |25j. The lesser 



Green's function was found from the Keldysh equation 
using the factor in the form (12). This means that 
is supposed to be unchanged due to the Coulomb 
interactions on the quantum dot, as it is proved to be 



within the Hartree-Fock approximation. We believe the 
proposed approximation to be appropriate at moderate 
and/or weak coupling of the dot to the leads when the 
coupling strengths are smaller than the characteristic 
Coulomb energy on the dot. We derived an explicit ex- 
pression for the occupation numbers of electrons in the 
dot (Eq. 12"!]) which enables us to compute them avoiding 
the long iterative procedure. The employed formalism 
gives correct results in the Coulomb blockade regime cor- 
roborating the results following from the transition rates 
(master) equations. Our approach is able to quantita- 
tively reproduce the relative heights (2 : 1) of the steps 
in the current voltage curves and to qualitatively describe 
the Kondo peak using the same expressions for the rel- 
evant Green's functions. To the best of our knowledge 
this was never achieved so far. 

Also, we remark that the present work addresses a gen- 
eral problem arising within the NEGF based computa- 
tions. To quantitatively describe fine features of Kondo 
effect and related phenomena one needs to include hy- 
bridization up to very high orders in calculations of the 
relevant Green's functions. Carrying out these cumber- 
some calculations where numerous approximations are 
inserted, it is easy to make a mistake and difficult to dis- 
cover one, especially in the absence of obvious small pa- 
rameters. Probably, such a mistake would not affect the 
very existence of the Kondo peak but it could distort its 
fine features. In the present work we propose a method 
which could help to verify the resulting expressions for 
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the Green's functions. We start from the clear point, 
namely, that the Green's functions suitable to describe 
the Kondo effect are suitable to describe the Coulomb 
blockade transport, as well. Therefore, one may verify 
obtained results by going to the Coulomb blockade limit. 
If it occurs that particular Green's functions fail to pro- 
vide the correct form of the I-V curves within this limit 
then one must conclude that there is some error in the 
adopted approximations for the Green's functions. 

It is shown that calculational scheme employed in the 
present work which uses a very simple approximation for 
G< brings quantitatively correct results for the electron 
transport through a quantum dot within the Coulomb 
blockade regime whereas advanced self-consistent cal- 
culations carried out in the recent paper [20| do not. 
This gives grounds to conjecture that it is not necessary 
(and, perhaps, it is not always correct) to use the same 
number of iterations in seeking approximations for re- 
tarded/advanced and the lesser Green's functions within 
EOM method. Also, we remark that validity of the 
Green's functions used in studies of the Kondo effect may 
be verified by applying them to calculate electron trans- 
port within the Coulomb blockade regime. The results 
could be generalized to include more realistic case of a 
dot including many electron levels, assuming that level 
separations are much greater than the Coulomb energy 
U, so one may neglect Coulomb interactions of electrons 
belonging to different levels. 
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APPENDIX 

Here, we explain how the expression (5) for the Green's 
function for the dot is derived. We introduce the fol- 
lowing notation for a general retarded Green's function 
Fourier component: 

((A, B)) = -i [ ({A(t)\B})e^ E+l, i^dt (26) 
Jo 

where A, B are operators, curly brackets denote the an- 
ticommutator and ( • ■ • ) stands for the average. For the 
Hamiltonian (1) the equation of motion for the retarded 
Green's function for the dot G r J(E) = ((c CT ; c+)) reads: 

(E — E + ij])((c a ;C+)) = 1 + ^ T r(j, a ({Crfia, C+}) 

r,/3 

+ U((n^c a -ct)). (27) 
Here, n_ CT = c^ a c- a , and the Green's function 
((crpcr', c t)) obeys the equation: 

(E - e r p a + iij)((c rl3a ;c+)) =T* pa ((c a ]C+)). (28) 

Substituting the expression for ((c r pa] c t)) determined 
by [2^1 into [271 we get: 

(E — Eq- E 0ct )((c ct ; c+» - 1 + f/((n_ ff c CT ; c+» (29) 

where the self-energy part £o<t is given by Eq. [5] 
The equation of motion for the four-operator function 
(in-tjC^c^^ includes higher order Green's functions. 
Omitting them we write out: 



{E-Eq-U + ir])((n- <7 c a ;4)) 

= {n-a) + { T rfj<r{{crl3an- a ; C+» + T r , 3 .- lT {{ct a C rl 3-- lT C lT ; C+» - T*p._ a { (c r& -- a C- a C a \ C+ ) )J . (30) 

r,0 

To proceed we must write equations for the Green's functions inserted in the right side of the Eq. [301 We get: 
(E - e rPa + if])((c rj3a n_ a ; c+)) 

^ T r/3a{{ n -<r C ^ C i)) + Yl fa Pi'* ( (c r&a cL a C r , & ._ a \ C+ ) ) - T* lp ._ a { fa P - a C^ a C r0VT ; C+) ) ) , (31) 
r'P 

(E - e r p.-„ + i77)((cl cr c r(3; _ .c .; c+)) 
= r *ft- < r(( rl -<rC -;c+» (r r ^_ .((cJ /3 ._ (7 c r/ 3 i ._ .c CT ;c+» +T r / (3CT ((ci CT c r(3; _ ff c r / (3a ;c+))), (32) 

r'P 

(E-2E -U + erft-v + irj)((c+p._ a c^c a ; c+» 
= - Tr/3;- a ((n-aC (T ;c+)) + ^ fa' p-,-* -*^^ c+ )) - r r / (3; _ (T ((c+ 3; _ o .c (T c r , / 3 ; _ (T ; c+») . (33) 

r'P 
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Writing out these equations [3T1I331 we neglected averages 
like ic^Crpa') including a creation/annihilation operator 
for the dot combined with the annihilation/creation op- 
erator for the source/drain reservoir. Such averages are 
omitted in further calculations, as well. Now, we decou- 
ple four-operator Green's functions included in the sums 
over "r"' in the equations following the scheme [2(| : 

((A+BC-D+)) « (A+B)((C;D+))-(A+C)((B;D+)). 

(34) 

Also, we use the approximation: 



(35) 



where f£ a is the Fermi distribution function correspond- 
ing to the energy e r/ 3 a . As a result, the Green's functions 
included in the right-hand side of the Eq. [30] get ex- 
pressed in terms of the Green's functions ((n_ CT c CT ; c+)) 
and ((ccr;c£:)). Substituting these expressions into Eq. 
[30] we obtain: 

Here, self-energy parts £20-, £q 2 are given by: 



r{3 



I Jr.— a 



l 



E 



17] 



E-2E Q -U + tr[)-a + ir\ 

(37) 



S 2<r = Y \ Tr ( } ~°\ 2 \ 
r0 



1 



E - 2E Q - U + e r/3 -„ + ir] 



1 



E - e r f3-a + if] 



^02 — ^Oo 



^2er- 



(38) 



(39) 



To arrive at the resulting expression for the dot Green's 
function G r J = ((c CT ; c+)) we substitute Eq. [35] into Eq. 
[Ml We have: 



G r J(E) 



(E-E - E 0a )(E -Eo-U- X£ 2 ) + UE la ' 

(40) 
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